Tunneling via individual electronic states in ferromagnetic nanoparticles 
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We measure electron tunneling via discrete energy levels in ferromagnetic cobalt particles less 
than 4 nm in diameter, using non-magnetic electrodes. Due to magnetic anisotropy, the energy 
of each tunneling resonance shifts as an applied magnetic field rotates the particle's magnetic mo- 
ment. We see both spin-increasing and decreasing tunneling transitions, but we do not observe the 
spin degeneracy at small magnetic fields seen previously in non-magnetic materials. The tunneling 
spectrum is denser than predicted for independent electrons, possibly due to spin- wave excitations. 

PACS numbers: 73.20.Dx, 75.60.-d, 73.23.Hk 



The forces that determine the electronic properties 
of semiconductor quantum dots or metal nanoparticles 
can be investigated in a particularly direct, fundamen- 
tal way through tunneling measurements of the discrete 
"electrons-in-a-box" spectrum of energy levels. This 
technique has been used to study the quantum-Hall-effect 
regime superconducting pairing in aluminum parti- 
cles and effects of more generic electron-electron 
interactions In this Letter we turn to a ferro- 
magnetic material (cobalt), with the aim of probing the 
ways in which strong exchange interactions and magnetic 
anisotropy affect the discrete electron spectrum, as well 
as investigating spin-polarized tunneling via single quan- 
tum states. Our work may be viewed as an extension to 
smaller size and/or lower temperature of previous exper- 
iments employing micron-size ferromagnetic islands ||] 
and nm-scale cobalt particles [||. We find a number of 
phenomena different from past studies of non-magnetic 
nanoparticles. We observe that the energy levels are cou- 
pled to the direction of the particle's total magnetic mo- 
ment, such that the levels shift as the moment is reori- 
ented. As might be expected, there is no degeneracy 
in Co between spin-up and spin-down tunneling levels 
near zero magnetic field, H. The energy spacing of the 
resonances is smaller than expected in an independent- 
electron model, suggesting the importance of low-energy 
many-body spin excitations. 

Our samples consist of Co particles connected to Al 
leads via tunnel barriers. The fabrication is similar 
to previous work on Al particles [|To|| , and a schematic 
sample cross-section is shown in Fig. 1(a). The top Al 
electrode is deposited first so as to fill a bowl etched 
through a SisN4 membrane (hole radius < 5 nm), and 
then a tunnel barrier is formed by oxidizing the Al. The 
Co nanoparticles are obtained by evaporation at room 
temperature of a Co layer 0.5 nm thick. Surface ten- 
sion causes the Co to form electrically separate particles. 
Scanning transmission electron microscope (STEM) im- 
ages of test samples in which 0.5 nm of Co is deposited on 
oxidized aluminum indicate Co particles with diameters 
1-4 nm, with center-to-center spacing 2-5 nm (Fig. 1(b)). 
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FIG. 1. (a) Cross-sectional device schematic, (b) Annular 
dark-field plan-view STEM image of Co particles. 

The crystal structure was not determined. The second 
tunnel barrier is formed on the Co particles either by de- 
positing a 0.8 nm-thick layer of Al at 77 K, which we 
oxidize in 50 mTorr of O2 for 3 min at room tempera- 
ture (sample 1), or by directly depositing 1.1 to 1.5 nm 
of AI2O3 (samples 2,3). Finally a thick layer of Al is 
deposited to make the second electrode. We select de- 
vices for which the cur rent- volt age curve at 4.2 K shows 
Coulomb staircase structure (not shown), indicating tun- 
neling via nanoparticles Jll]] . 

Figure 2 shows the tunneling spectra at the onset of 
conduction for the first Coulomb threshold, for three 
samples. The spectra consist of well-resolved peaks due 
to tunneling via discrete electronic levels within each 
particle, qualitatively similar to previous measurements 
in Al and Au 0,0] . The values on the energy axis are de- 
termined by dividing the voltage by (C1+C2VC2, to cor- 
rect for capacitive division of the bias, where Ci and C2 
are the capacitances of the particle to the two electrodes 
[ p!o| . We can determine this ratio to within 20% by fitting 
the temperature-dependent broadening of peaks ]T^| , or 
in one sample having no voltage-dependent charge shifts 
(sample 2), we can achieve 1% accuracy by comparing the 
voltages required for tunneling via the same electronic 
states at positive and negative bias The peak spac- 
ing for all three samples is much less than the Coulomb 
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FIG. 2. Tunneling spectra of 3 different samples at T=20 
mK and fi^H = 1 T . H is parallel to the SiaN4 membrane. 
The energy is obtained by dividing V by (C1+C2VC2 = 1.2, 
2.17 and 2.5, respectively, for the three samples. 

charging energy (> 30 meV, based on particle size), so 
that all peaks in each spectrum correspond to tunneling 
via states with the same number of electrons, either one 
more or one less than the initial state. 

Unlike non-magnetic particles where the energy levels 
have a simple linear dependence on H due to the spin Zee- 
man energy [ p0[ , the levels in Co have a strong, non-linear 
dependence on H for small H. Figure 3(a) shows the en- 
ergy of the first three tunneling resonances of sample 1, as 
jhqH was swept about a hysteresis loop from —0.45 T to 
0.45 T and back again. Starting at —0.45 T (thick lines), 
the tunneling energies shift in a continuous manner as the 
field is ramped to zero and the magnetic moment vector 
m of the particle relaxes toward its easy direction. As the 
field is ramped further, we observe a sudden jump in all 
three transition energies, at fioH sw = 0.23 T. We inter- 
pret this jump as due to the reversal of m in the single- 
domain Co nanoparticle [|l3[ . The energy shifts with H 
are hysteretic with respect to the direction of the field 
sweep, with the expected symmetry around H = 0, and 
the scans are repeat able over days. The value of H sw is 
comparable to SQUID measurements of 25 ± 5 nm diam- 
eter Co particles JyJ, and corresponds to an anisotropy 
energy density of order K=10 5 J/m 3 . 

These curves indicate a significant coupling between 
the level energies and the orientation of the magnetic 
moment of the nanoparticle. Let us consider the sim- 
plest model for this anisotropy, so as to see what fea- 
tures of the data may be explained simply, and what 
features may require deeper understanding. We call 
the operator for the total electronic spin with N elec- 
trons HS(N), and we assume that the anisotropy and 
Zeeman energies are sufficiently weak relative to the 
exchange splitting between different spin multiplets 
that the magnitude S in the ground state remains 



0.6 -1 — 1 1 1 1 1 1 1 1 r 




FIG. 3. (a) Hysteresis curves showing the dependence of 
tunneling energies on H for sample 1, at T=20 mK. (b,c) 
Lowest-energy transitions calculated using the Hamiltonian 
discussed in text, for S=50 and H oriented 45° from the easy 
axis, for the case where S increases during tunneling (b) and 
decreases (c). The qualitative features are independent of the 
value of S. (d) Voltage threshold for tunneling in one sample 
which exhibits anti-hysteretic behavior. 

constant as H is varied, so that we can perform our cal- 
culation in the space of this one spin multiplet. (In a 
classical language this corresponds to assuming that the 
ground state magnetic moment fh(N) simply rotates as 
a function of H.) Including the Zeeman energy and the 
simplest model of easy-axis anisotropy in the z direction 
]l5[ , the i^-dependent terms in the Hamiltonian within 
the ground-state multiplet for N electrons can be written 

H = -g^B^H • S - K m ii B S 2 z />JS{S + l) (1) 

where K m is an anisotropy energy per unit \m\. 
(K rn =K ( volume)/ 1 m\.) We have omitted from Eq. 1 the 
charging energy, which we assume to be independent of 
H. For given total spin values in the N and 7V±l-electron 
states, we diagonalize the Hamiltonian numerically to 
find the energy levels, and then calculate the allowed tun- 
neling transition energies as E(N±1 , H) — E(N , H). We 
allow the moment vector of the nanoparticle to undergo 
reversal at the classical switching field, which depends 
on the angle between H and the easy axis. Represen- 
tative results for the lowest-energy tunneling transitions 
are shown in Fig. 3(b) for the case that S increases dur- 
ing tunneling, and in Fig. 3(c) for S decreasing. The 
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model successfully reproduces the hysteretic energy max- 
ima near zero field observed in the lowest energy transi- 
tion (with S increasing), as well as the existence and the 
sign of the abrupt switching to lower energy at H sw . If 
we define Sh as the component of the total electron spin 
in the direction of H, in general we find that tunneling 
transitions in which \(Sh)\ increases give maxima near 
H = and transitions decreasing \(Sh)\ give minima. We 
can identify both types of behavior in Fig. 3(a). We have 
also solved the classical analogue of the model, which 
gives similar results for the ground-state to ground-state 
transitions. 

The measured tunneling energies often have small-scale 
non-monotonic variations as a function of H that are not 
present in our minimal model. A likely cause is magnetic 
interactions between nearby Co particles. The most dra- 
matic example we have observed is shown in Fig. 3(d), 
where we plot the independence of the threshold volt- 
age for tunneling via a Co nanoparticle too large for dis- 
crete resonances to be observed. We see anti-hysteresis 
- magnetization reversal occurs before H changes sign. 
This can be explained by the influence of a dipolar mag- 
netic field oriented opposite to the applied ii", produced 
by a second magnetic nanoparticle adjacent to the one 
through which electron tunneling occurs. The reversed 
field from the second particle can shift the hysteresis 
curve of the first so that its value of H sw can be neg- 
ative, while the non-monotonic shifts at large positive H 
(0.2 T) are understood as the magnetization reversal of 
the second particle. A dipole field 5 nm from the cen- 
ter of a 3-nm-diameter Co particle is of order 0.1 T, so 
that Zeeman interactions alone are not strong enough to 
explain the level shifts that we observe. 

A second departure of the data from the minimal 
model is that in some cases the model fails to describe 
the combined low and high H variations of the transition 
energies (Fig. 4). At large ii", the simple model predicts 
a linear extension of the low field curves, with the tunnel- 
ing energies Zeeman-shifting as ±g e fi/j,BlJ>oH/2. Figure 
4(a) shows that at large H the ground state transition 
of sample 1 moves to higher energies with increasing 
H, indicative of a tunneling transition in which \(Sh)\ 
decreases, whereas the low-field behavior indicates a 
\(Sh)\ -increasing transition. Similarly, the ground state 
transition of sample 2 (Fig 4(b)) is \(Sh)\ -decreasing 
at low field (dip to lower energy) but \(Sh)\ -increasing 
at high field (energy shift to lower energies). We see 
two possible explanations. Either (i) the true form of 
the anisotropy is more complicated than assumed in 
Eq. (1), such that K m fluctuates to have different val- 
ues for different spin multiplets, or (ii) contrary to our 
model's other initial assumption, the total spins of the 
ground-states for N and N ± 1 electrons are not in- 
dependent of H, so that S (classically, m) does not 
simply rotate as H is ramped. The second possibility 
would mean that the spin character of the many-body 




FIG. 4. Tunneling energies over a larger range of H. 

electron ground states could evolve with ii", so it is pos- 
sible that the threshold events for tunneling may consist 
of \(Sh)\ -increasing transitions for some values of H, and 
\(Sh)\ -decreasing transitions elsewhere. 

For /iqH > 2 T, all measurable transition energies 
within a given sample have the same sign of shift as a 
function of H. This is different from results in Al and 
Au nanoparticles, for which Zeeman spin splitting of each 
orbital state gives rise to both upward and downward- 
trending states vs. ii", with comparable conductance am- 
plitudes, and with a degeneracy at ii" = 0. The ab- 
sence of spin degeneracy is not surprising in Co, since 
the exchange field breaks the symmetry between spin-up 
and down. However, since Co is not fully spin-polarized 
(P w 30%, [fl6||), the absence of observed peaks shift- 
ing in both directions is an indication that electrons do 
not couple equally well to all allowed many-body energy 
levels. One factor may be that Clebsch-Gordan coef- 
ficients will be very different for \(Sh)\ -increasing and 
\(Sh)\ -decreasing transitions in a system with a large 
ground-state value of the total electron spin, so that some 
tunneling matrix elements may be immeasurably weak. 
Note that there is a sign difference between the large- ii" 
slopes for sample 2 and those of samples 1 and 3. This 
might reflect a difference between electron-addition and 
electron-removal tunneling processes, a possibility which 
could be tested by fabrication of gated devices. 

We wish to distinguish the linear dependence on ii" 
that we measure for the tunneling resonances at large H 
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from the linear shift in the chemical potential of micron- 
scale Co islands described by Ono et al. ||. The data of 
H can be understood in terms of level crossings in a mag- 
netic island with a continuum of energy states. As H is 
ramped, spin- up levels cross with spin-down on account 
of the Zeeman energy, and because the densities of the 
spin levels are different, more states move one way than 
the other and the chemical potential must shift. This 
model is not applicable in explaining the shifts of the in- 
dividual levels that we observe, because in our case any 
level crossings should be individually resolvable. Quan- 
titatively, the shifts examined by Ono et al. in Co corre- 
spond to an effective g-factor of 0.7 ||, smaller than the 
values we measure. (See Fig. 4.) 

Finally, we turn to the measured density of tunneling 
resonances shown in Fig. 2. For the range of particle 
sizes imaged by STEM, 1-4 nm in diameter, the average 
level spacing predicted for simple non-interacting elec- 
trons should be between 0.75 and 40 meV, given the 
calculated density of states (including both sp and d 
bands) in Co of 0.88 eV^atom- 1 0. In 

our measured 

spectra, the energy spacing between tunneling peaks is 
less than 0.2 meV. Enhanced densities of tunneling reso- 
nances have previously been seen in Al nanoparticles for 
values of voltage much greater than the single-electron 
level spacing, and they were explained as an effect of non- 
equilibrium electron-hole excitations within the particle 
||. The data in Fig. 2 for samples 1 and 2 are differ- 
ent from the Al results, however, in that an increased 
density of levels is observed for energies even below the 
expected single-electron spacing, where electron-hole ex- 
citations should not be produced during tunneling. As 
explanation, we note that electron excitations within a 
Co nanoparticle may include low-energy spin waves in ad- 
dition to the independent-electron-type excitations seen 
in Al. Inelastic emission of spin waves during tunnel- 
ing may directly contribute new tunneling peaks, and/or 
non-equilibrium spin excitations generated during tun- 
neling may produce extra tunneling peaks by shifting 
single-electron states As a check, we can estimate the 
minimum energy needed to excite spin- wave modes. For 
a spatially uniform mode (fc=0), the excitation energy 
can be calculated using the anisotropy term in Eq. (1), 
to be « 2K m \iB' If we use the size of the jump in the 
tunneling energy at H sw , AE w 0.05 meV (Fig. 3) to es- 
timate K m (using AE « /i#i<r m , see Fig. 3(b,c)), we ar- 
rive at a value of 0.1 meV for the spin- wave energy, which 
would explain the enhanced density of tunneling states. 
The contribution of exchange energy to the lowest-energy 
non-uniform spin- wave modes can be estimated by quan- 
tizing the spin- wave dispersion curve of Co within the size 
of a nanoparticle. This gives an energy (300 meV)(a/d) 2 
where a is a lattice spacing and d is the particle diameter 
[|l8f , or « 1 meV for a 4-nm particle. 
In conclusion, we have measured discrete tunneling res- 



onances in nm-scale ferromagnetic Co particles. Mag- 
netic anisotropy causes each resonance energy to shift 
reproducibly by on the order of 0.1 meV as H is swept 
about the hysteresis loop. This effect may provide a 
means to probe the dynamics of magnetization reversal in 
nm-scale particles, complementary to magnetic force mi- 
croscop y ]l9| ], Hall magnetometry ]20| , and SQUID tech- 
niques [[I4|. Qualitative features of these shifts can be 
described by a simple model. However, a full explanation 
of the measurements will require a more detailed under- 
standing of the electronic states inside a ferromagnet, in- 
cluding at least the contributions of low-energy collective 
spin- wave excitations to the electron states, and effects 
of Clebsch-Gordan coefficients in tunneling. 
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